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ABSTRACT 

Spatially discrete difference approximations of hyperbolic initial-boundary- 
value problems (IB VP’s) require numerical boundary conditions in addition 
to the analytical boundary conditions specified for the differential equations. 
Improper treatment of a numerical boundary condition can cause instability 
of the discrete IBVP even though the approximation is stable for the pure 
initial-value or Cauchy problem. In the discrete IBVP stability literature 
there exists a small class of discrete approximations which are called border- 
line cases. For nondissipative approximations, borderline cases are unstable 
according to the theory of the Gustafsson, Kreiss, and Sundstrom (GKS) but 
they may be Lax-Richtmyer stable or unstable in the L 2 norm on a finite 
domain. We show that borderline approximations can be characterized by 
the presence of a stationary mode for the finite-domain problem. A station- 
ary mode has the property that it does not decay with time and a nontrivial 
stationary mode leads to algebraic growth of the solution norm with mesh 
refinement. We give an analytical condition which makes it easy to detect a 
stationary mode, and we investigate several examples of numerical boundary 
conditions corresponding to borderline cases. 

1. INTRODUCTION 


In this paper we consider the stability of spatially discrete approximations to hyper- 
bolic IBVP’s. For simplicity we consider the stability of approximations to the IBVP 
for the model hyperbolic equation 


du _ du 
dt C dx' 


0 < x < L, 


t > 0 


( 1 . 1 ) 


where c is a real constant. One must specify initial data at t = 0, and the IBVP is 
well-posed if an analytical boundary condition is prescribed at x = L 


u(L,t) = g{t) for c > 0. , (1.2) 

A semidiscrete approximation of (1.1) is obtained by dividing the spatial interval 
into J subintervals of length Ax where J Ax = L, x = Xj = jAx and approximating 


1 Reprint for the Second International Conference on Hyperbolic Problems, Aachen, Fed. Rep. of 
Germany, March 14-18, 1988 



the spatial derivative u x by a difference quotient. As a prototype approximation, we 
replace u z by a second-order-accurate central-difference quotient to obtain a system of 
ordinary differential equations (ODE’s) 

= 2 Kx ^ Uj+1 ~ Uj_1 )’ J = 1 ’ 2, ' ,J_1 ( L3 ) 

where u ; = uj(t) denotes the approximation to tx(x, t). The right boundary (x = L) is 
advanced by using the analytical boundary condition (1.2). For the stability analysis, 
we assume that the boundary condition is homogeneous, i.e., g(t) = 0, and for the 
semidiscrete problem we write uj = 0. 

A complication in completing the approximation is the fact that more boundary 
conditions are required for the semidiscrete approximation than are specified for the 
partial differential equation. If we apply (1.3) at the left boundary (j = 0), then the 
difference stencil protrudes one point to the left of the boundary. It is clear that that 
a numerical boundary condition is required. For example, at the left boundary (j = 0) 
we can change from a centered approximation to a one-sided approximation of u x : 

^ 7 - = ~ 7 -[— Ortt 2 + (1 + 2a)uj - (1 + <*)u 0 ] (1.4) 

at Ax 

where a is a parameter. Any procedure, e.g., (1.4), used to provide a numerical bound- 
ary condition is called a numerical boundary scheme (NBS). In the stability analysis of 
this paper it is convenient to express the NBS as a space extrapolation formula. The 
NBS (1.4) is equivalent to 


g(£)u_! = 0, where q{E) = (E - 1) 2 (2 aE - 1) (1.5) 

and the shift operator E is defined by 

E u j — Uj . (l*6) 

The system of ODEs (1.3) together with the analytical boundary condition uj = 0 
and the NBS (1.4) can be written in vector-matrix form as 

~ = -MO (1.7) 

where u is a ./-component vector and A is a J x J matrix. The essential element in 
the stability of the semidiscrete approximation (1.7) is the behavior of the solution at 
a fixed time as the spatial mesh is refined. Consequently, one must consider an infinite 
sequence of ODE systems of dimension J where J — ► oo as Ax — > 0. 

For the semidiscrete approximation (1.3) with periodic boundary conditions, the 
matrix A = A p is a skew-symmetric circulant matrix. Consequently the eigenvalues of 
A p are pure imaginary and the semidiscrete approximation is said to be nondissipative. 
The analysis of stationary modes for dissipative approximations will be given in a 
subsequent paper. 
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Stability of a semidiscrete approximation with homogeneous boundary data means 
that there exists an estimate of the solution in terms of the initial data. For example, 
the semidiscrete approximation represented by the sequence of ODE’s (1.7) is Lax- 
Richtmyer stable if there exists a constant K > 0 such that for any initial condition 
u(0) 

NON < *llu(0)l| (1.8) 

for all J > 0, J Ax = L and for all t, 0 < t < T with T fixed. In this paper the symbol 
|| • || denotes the discrete Li norm defined by 


I|u(0ll2 



1/2 


(1.9) 


Two methods for carrying out a stability analysis are the energy method and the 
normal mode analysis. The normal mode analysis is an eigenvalue analysis. If we look 
for a solution of (1.7) of the form u(<) = then we obtain Aj — s$. But this is 
just the eigenvalue-eigenvector problem for the matrix A where ^ is the eigenvector 
and s is the eigenvalue. The practical problem of implementing stability tests on the 
eigenvalues is that the normal mode analysis for a discrete hyperbolic IBVP on a finite 
domain is, in general, analytically intractable. 

The intractability can easily be demonstrated by the normal mode analysis of the 
ODE system (1.7). The components <j>j of the eigenvector <f> and the normalized eigen- 
value s = (A x/c)s are given by 

4>j = a\t c J - (— k?) j (— l/>c) J ], 2s = k — ~ (I.10a,b) 

K 

where a is an arbitrary nonzero constant and k is a root of the characteristic equation 

q(ic) - (— K. 2 ) J + 1 q{— 1/k) = 0. (1.10c) 

The polynomial g(/c) depends solely on the NBS, i.e., q(K.) is the polynomial associated 
with the NBS written as an extrapolation formula. For example, the polynomial q(n) 
for the NBS (1.5) is 

q(K.) — (k — l) 2 (2a/c - 1). (1.11) 

Since J Ax = L, the degree of (1.10c) increases as the spatial mesh increment Ax 
decreases. In general, one cannot solve for the roots of (1.10c) which accounts for the 
analytic intractability of the normal mode analysis on a finite domain. 

Although the Lax-Richtmyer condition (1.8) is a conventional stability definition, 
there is no known algebraic test to check the stability condition for discrete hyperbolic 
IBVP’s on a finite domain. A more complicated stability definition is used in the 
theory developed by Gustafsson, Kreiss, and Sundstrom (GKS) jl). Strikwerda [4] 
has extended the GKS theory to semidiscrete approximations. The advantage of the 
GKS theory accrues from the fact that a finite-domain problem with two boundaries 
is divided into a Cauchy problem and two quarter-plane problems each of which can 
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be analyzed separately by a normal mode analysis. The analogues of (1.10) for the 
right-quarter plane problem are 




(I.12a,b) 


where k is a root of 

q(n) = 0. (1.12c) 

This is the same polynomial gr(/c) that appears in (1.10c). The roots of (1.12c) are 
easily found since q(t c) is of low degree. Algebraic tests based on the roots of q(n) and 
the corresponding eigenvalues s provide necessary and sufficient conditions for GKS 
stability. 

The connection between the normal mode analysis for the finite-domain problem 
and the normal mode analysis for the quarter-plane problem is rather obscure. In 
a recent paper [6], we derived asymptotic estimates of the eigenvalues for the finite- 
domain problem. These estimates were used to relate the normal mode analysis of 
the finite-domain problem and the GKS quarter-plane analysis. In order to derive the 
asymptotic estimates for the roots of (1.10c), we assumed that particular roots can be 
identified for each J and we write k = k(J). Furthermore, there is no loss of generality 
in assuming that \k.(J)\ < 1. We showed that the roots of the characteristic equation 
(1.10c) can be divided into three distinct classes according to the asymptotic behavior 
of \k(J)\ j in the limit J —* oo. For |»c(J)| < 1 there are only two possibilites: 

(I): lim |/c(J)|' / = 0, (II) : lim |K(J)j J = constant > 0. (1.13) 


For roots in class (/), it is clear that (1.10c) reduces to the quarter-plane equation 
(1.12c) as J — * oo. Consequently, a root in class (I) becomes a root of the quarter- 
plane polynomial (1.12c) in the limit J — > oo. The eigenvalues corresponding to the 
/c’s of class (II) are benign in the stability analysis in the sense that they satisfy the 
inequality 


»(s) < 


const 


J 


00 


(1.14) 


which is a necessary condition for stability (see [6]). 


2. STATIONARY MODES 


In addition to (I) and (//), there is a third class of roots 


(///) : |ic(J)| = 1. 


( 2 . 1 ) 


If /c is in class (III) and is independent of J, then k. remains fixed on the unit circle 
for all J. For this to happen it is obvious from (1.10c) that the polynomials q(ic) and 
q(— l//e) must have a common factor. This common factor leads to identical roots for 
both the quarter plane polynomial g(/c) and the finite-domain characteristic equation 
(1.10c). These roots are fixed (independent of J) on the unit circle and from ( 1.10b) 
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one obtains 9?(s) = 0 with S(s) = fixed. Consequently, there is a stationary mode, i.e., 
a mode with k and s independent of J. 

If there is a stationary mode for the finite-domain problem, the GKS perturbation 
test will always indicate the presence of a GKS generalized eigenvalue and consequently 
the semidiscrete approximation is GKS unstable. The existence of a GKS generalized 
eigenvalue is easily proved. Since |<c| = 1 and both k and — 1/k are roots of q(n,) = 0, 
we need two perturbation tests, i.e., 

K a = I 1 ~ £ ) e ’*> K b = “(I - «)«“** (2-2) 

where /c* and /c£ denote the perturbations of k and —1/k inside the unit circle. But 
from (1.10b) it follows that 


*{*[<)) = -»[■§«)] # 0 (2.3) 

and consequently there is a generalized eigenvalue 3?(i *)> 0 corresponding to either k 
or —1/k. 

The importance of a stationary mode is the following. Gustafsson et al.[2] have 
proved that if the Cauchy stability requirement of the GKS theory is replaced by a 
more stringent energy estimate, then GKS stability implies Lax-Richtmyer stability in 
the L 2 norm. There is a small number of known examples showing that Lax-Richtmyer 
stability in the L 2 norm does not imply GKS stability. These examples all involve what 
are called borderline cases. We show that borderline cases have a stationary mode for 
the finite-domain problem. The GKS quarter plane analysis cannot detect whether or 
not a particular mode is stationary for the finite-domain problem. However, from our 
analysis, stationary modes are easy to detect since they occur if and only if g(«) and 
g(— l//c) have a common factor. 

From the point of view of an eigenvalue analysis, a semidiscrete approximation with 
a stationary mode must be treated separately since any instability derives not from an 
eigenvalue with a positive real part but from the algebraic growth (as J — > 00 ) of the 
norm of the solution. 

3. EXAMPLES OF STATIONARY MODES 

In this section we examine several examples of stationary modes arising from various 
NBSs for semidiscrete approximations. We analyze the first example in detail but give 
only a summary of the stationary mode analysis for the subsequent examples. All of 
the examples are GKS unstable. Examples 1 and 2 are Lax-Richtmyer unstable but 
examples 3 and 4 are Lax-Richtmyer stable. In each of the examples we follow the 
convention of having the boundary condition of interest on the left boundary, i.e., the 
GKS analysis is done for a right quarter-plane problem. 

3.1 EXAMPLE 1 

Our first example of a semidiscrete approximation with a stationary mode is NBS 
(1.4) or equivalently (1.5) for a = —1/2. The polynomial (l.ll) becomes 

?(k) = -(/c 2 - l)(/c - 1), </(— l//c) = — (*c 2 — 1 )(k + 1)k 3 . (3.1) 
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The polynomials g(/c) and q(— 1 /k) have the common factor ( k 2 - 1), and consequently, 
there is a stationary mode. The characteristic equation (1.10c) has the roots 


K, — il 


(3.2) 


independent of J and, from (1.10b), s = 0. According to a GKS stability analysis, a = 
— 1/2 is the borderline case between stability (a > —1/2) and instability (a < -1/2). 
For a = —1/2 there is a stationary mode, and consequently, this borderline case is 
GKS unstable. 

The NBS (1.4) with a = -1/2 is 


du 0 

~dt 


c 

2Az 


(u 2 - Uo). 


(3.3) 


As an aid to interpreting the solution of a semidiscrete approximation with (3.3) as the 
NBS, it is useful to consider zeroth-order extrapolation 


Uo{t) = Ui(t) 


(3.4) 


as an NBS. Differentiation of (3.4) with respect to time yields 

duo du\ 

dt dt 


(3.5) 


By replacing du\jdt in the above equation by the interior approximation (1.3) evaluated 
at j = 1, one obtains the NBS (3.3). If we integrate (3.5) from 0 to t then 


u 0 (t ) — ui(t) = u o (0) — Uj(0) = constant. (3.6) 

Hence, the NBS (3.3) is equivalent to zeroth-order space extrapolation if the initial 
value is reset at the outflow boundary so that 


ti 0 (0) - «i(0) = 0. 


(3.7) 


A semidiscrete approximation with zeroth-order extrapolation (3.4) as the NBS 4s 
both GKS stable and Lax-Richtmyer stable. However, differentiation of the NBS (3.4) 
yields an approximation which is both GKS unstable and Lax-Richmyer unstable. Next, 
we show that the Lax-Richtmyer instability is due to the presence of a stationary mode. 

The stationary eigenvector corresponding to the stationary roots (3.2) is found by 
substitution of either « = 1 or — 1 into (1.10a) to obtain 


<t>j = a[l-(-l)* +J ], j = 0, 1, , • • • , «7 — 1 (3.8) 

where a is an arbitrary nonzero constant. The L 2 norm of the eigenvector (3.8) is 


\\*h 



\a\y/2L, 

|a|v/2L(J + l)/J, 


J even 
Jodd 


(3.9) 
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where we have used J Ax = L. 

The matrix A corresponding to the ODE system (1.3) with (3.3) as the NBS and 
u j = 0 as the analytical boundary condition is A = (c/ Ax)A where 




o 

o 


o 


o 


1 0 1 

- 1 0 


(3.10) 


Since the first two rows are equal, it is obvious that the matrix A has an eigenvalue s = 
0. The eigenvector (3.8) is the right eigenvector ^ of the matrix (3.10) corresponding 
to s = 0. The left eigenvector £ corresponding to s = 0 is easily found by inspection of 
the matrix (3.10) to be 

e r = (-1,1,0, - ,o]. ( 3 .H) 

If we choose a = 1/2 in (3.8), then the inner product of £ and ^ is 

* r *=±l (3.12) 


where the plus sign is used if J is even and the minus sign if J is odd. 

The matrix (3.10) has a complete set of eigenvectors and consequently the general 
solution of the ODE system (1.7) can be written as 


u(t) = V 0 t e-m 


£—0 


(3.13) 


where ^(£) denotes the eigenvector corresponding to the eigenvalue sg. We denote the 
stationary eigenvector (3.8) by ^(0) and rewrite (3.13) as 

J - i 

u(o = m(o) + y. w'm- (3-i4) 

£=1 

The sum on the right hand side of (3.14) consists of the eigenvalues and eigenvectors 
for zeroth-order space extrapolation (3.4) as the NBS. 

For given initial data, the coefficient f3 0 associated with the stationary eigenvector 
^(0) is determined by taking the inner product of the left eigenvector £ given by (3.11) 
with u(0) given by (3.14): 


£ r u(0) = [-uo(0) + iMO)] = ±p 0 (3.15) 

where the plus sign corresponds to J even and the minus sign to J odd. In obtaining 
(3.15) we have made use of the fact that the right and left eigenvectors of a matrix are 
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orthogonal. Formula (3.15) for the coefficient fio is consistent with our earlier assertion 
that the NBS (3.3) is equivalent to zeroth-order space extrapolation if the initial value 
is reset at the outfow boundary so that Uo(0) — ui(0) = 0. 

One can show analytically that 5R(s£) < 0 for all the modes except the stationary 
mode, and consequently, the asymptotic solution ( t —* oo) is from (3.14) 


u (0 « 0), 

The 1/2 norm of (3.16) for J even is 


oo. 


M<)|| 2 « M*(o)||, « l/fcly'l = K(0) - «,(o)|^|. 


OO 


(3.16) 


(3.17) 


where we have used (3.9) with a = 1/2 and (3.15). 

The Lax-Richtmyer instability resulting from the NBS (3.3) arises from the presence 
of the stationary mode. We illustrate this on the finite domain 0 < x < L by choosing 
the initial data 

1 for 7 = 0 

, (3.18) 

0 for j > 0 


«>( 0 ) = 

where J Ax = L. The Li norm of this initial data is 

mo)ii 2 = a* 1/2 = 


(3.19) 


For the initial data (3.18) with J even, the Li norm of the asymptotic solution (3.17) 
is 

HOII2 


From (3.20) and (3.19) it follows 

NOIIa 


VJ 

V2 
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T 

(3.20) 

||u(0)|| 2 , t -4 OO. 

(3.21) 


Consequently, the Li norm of the solution is not uniformly bounded on 0 < t < T (with 
T large) for the initial data (3.18) and the semidiscrete approximation is Lax-Richtmyer 
unstable. 

If the initial data are smooth, then the semidiscrete approximation converges al- 
though the global order of accuracy drops to first order. If we assume that the initial 
data are smooth , then 


«o(0) - uj(0) = /(0) - /(Ax) = -Ai— ^ + 0(Ax 2 ). 


In this case (3.17) becomes 

IMOIk 


0/(0) 


dx 


— Ax = O(Ax), t -4 00 . 


(3.22) 


(3.23) 
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3.2a EXAMPLE 2a 


In this section the model hyperbolic equation is 
du du 

— = -c — , c > 0, 0 < x < L, t > 0, (3 24) 

dt dx ’ 

i.e., the analytical BC is specified at the left boundary, x = 0. The spatial derivative 
u x is approximated by a centered approximation and the PDE (3.24) is replaced by 
the system of ODEs 

du, c , 

*- = 2Ai (Uj - 1 " U) ' +,) ’ ( 3 - 25 > 

The homogeneous analytical boundary condition is 


u 0 = 0. 


(3.26) 


The NBS is 

duj c , , 

IT = Tr (Uj -' ‘ Uj} - < 3 ' 27 > 

The semidiscrete approximation (3.25) with boundary conditions (3.26) and (3.27) 
is both Lax-Richtmyer and GKS stable. In example 1 we showed that differentiation 
with respect to time of a stable NBS, i.e., zeroth-order extrapolation, resulted in an 
unstable semidiscrete approximation. In this section we show that differentiation of 
the analytical boundary condition (3.26) with respect to time, i.e., 


du 0 

dt 


= 0 


(3.28) 


leads to an unstable approximation. The source of the instability is the introduction of 
a stationary mode. It should be noted that this example differs from the other examples 
in this paper since the stationary mode is introduced by an improper modification of 
the analytical boundary condition, i.e., replacement of (3.26) by (3.28), rather than an 
improper choice of the NBS. 

It is easy to show by evaluating (3.25) at j = 0, that (3.28) is equivalent to the 
extrapolation formula 


h(E)u - j = 0, where h(E) = 1 — E 2 . (3.29) 

Since Ii(k) = 1 — /c 2 it follows directly that 

h(- 1/k) = — /i(k) / k 2 (3.30) 

as one can easily verify. Therefore, /i(k) and h(- 1 /k) have a common factor, namely 
h(t c) itself, and consequently there is a stationary mode. The roots of fi(#c) = 0 are /c = 
±1. The stationary mode has the eigenvalue s = 0 and the corresponding eigenvector 
is 

<t>j = a, j = 0,1,2, • -,J (3.31) 
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i.e., a vector with constant elements. 

Since all the modes except the stationary mode are damped for large time t, the 
asymptotic solution is 

u(t) « t • oo (3.32) 

where /? 0 = uo(0) and the elements of ^(0) are given by (3.31). For the initial data 
(3.18) the asymptotic solution is 

HOIh « \fL = \/j||u(0)|| 2 t -► oo. (3.33) 

Consequently, the norm of the solution is not uniformly bounded for 0 < t < T (with T 
large) for the initial data (3.18) and the semidiscrete approximation is Lax-Richtmyer 
unstable. 

3.2b EXAMPLE 2b 

We return to the model equation (l.l) with the NBS on the left boundary, x = 0. If 
we choose the overspecified Dirichlet condition 

u 0 = 0 (3.34) 


as the NBS, the resulting semidiscrete approximation is both Lax-Richtmyer and GKS 
stable. However, if we differentiate (3.34) with respect to time 


du o 
dt 


= 0 


(3.35) 


the resulting approximation is both Lax-Richtmyer unstable and GKS unstable. 

The NBS (3.35) is equivalent to the extrapolation formula 

q(E)u_ i = 0, where q(E) = 1 - E 2 . (3.36) 


The polynomial q(n :) = 1 — #c 2 , and therefore, 

g{- 1 /k) = -?(k)/k 2 - 


(3.37) 


Consequently, there is a stationary mode. For this example and the following examples, 
3 and 4, all of the roots of the characteristic equation are in class {III), i.e., \k{J)\ = 1 
and all the eigenvalues are pure imaginary. Consequently, the normal modes for the 
finite-domain problem can be found analytically. For J even the stationary mode has 
a repeated eigenvalue s = 0. For sufficiently large time the asymptotic solution is 

( — j — odd 

u i(0 w ^ L , t — ► oo. (3.38) 

(0, j = even 

For the initial data (3.18), the norm of the asymptotic solution is 

HOII 2 « ^^H°)ll2, «- 00. (3.39) 
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The Lax-Richtmyer instability arises from the factor \fj and not from the linear growth 
in t. 


3.3 EXAMPLE 3 


Consider the inconsistent NBS 


u_ i - — Ui 


(3.40) 


which can be written as 


q(E)u-i = 0, where q(E) = 1 + E 2 . (3-41) 

This NBS was used by Trefethen [5] for the leap frog scheme as an example of a fully 
discrete approximation to (1.1) that is Lax-Richtmyer stable and GKS unstable. As 
outlined below the same stability results are obtained when the NBS (3.40) is used for 
the semidiscrete approximation (1.3). 

Since g(/c) = 1 + /c 2 , one has 


?(-!/*) = 9(«0/k 2 (3.42) 

and consequently, there is a stationary mode. From ( 1.10c) one obtains 

q( K )[l + (~k 2 ) j ] = 0. (3.43) 

The roots of the equation (3.43) are determined from 

q(K.) = 1 + k 2 — 0, and 1 + (-/c 2 ) J = 0. (3.44a, b) 

For the stationary mode the roots of q(t c) = 0 are k = ±* and the corresponding 
eigenvalues are s = ±t. However from (1.10b) the roots k = ±t both lead to trivial 
stationary eigenvectors for the finite-domain problem. 

For this example the eigensolutions of the finite-domain problem can be found ana- 
lytically. The remaining roots of the characteristic equation can be determined from 
(3.44b) by using the roots of unity formula. The corresponding eigenvalues are pure 
imaginary. In particular, for J odd the eigenvalues are 

/ 2 ^ V 

si = isiny—j-J, £ = 1,2,--,J (J odd). (3.45) 

Here we have the rather amazing result that for J odd, the eigenvalues are analytically 
identical to those of the circulant matrix associated with the spatially periodic problem. 
The eigenvectors are, of course, different. 

One can show by the energy method that the semidiscrete approximation is Lax- 
Richtmyer stable. In particular 

MOIh < V5 ||u(0)|| 2 . (3.46) 
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On the other hand, the approximation is GKS unstable because there is a stationary 
mode. 

We briefly outline the GKS perturbation analysis. For k — i, the eigenvalue is s = i 
where obviously 9t(s) = 0 and hence we must check to see if there is a GKS generalized 
eigensolution. Let k* denote a perturbation of k = i which is inside the unit circle, i.e., 

k = e ie (l - e), #=- + €, £ > 0. (3.47) 

2 


By inserting (3.47) into (1.10b), we find 

SR(s*) = e 2 + 0(e 3 ). (3.48) 

Since the perturbation of 2R(J) is positive there is a GKS generalized eigenvalue and 
semidiscrete approximation is GKS unstable. The fact that an e perturbation of k 
yields an f 2 perturbation in 9J(s) indicates a weaker type of GKS instability than for a 
conventional GKS generalized eigenvalue where the perturbation in k and $t(s) are of 
the same order. 


3.4 EXAMPLE 4 

Our last example is due to Gustafsson[3]. The wave equation, Ut< = u x *, written as 
a first-order system is 



u 


— i 

o 

t— ‘ 

w = 


, B = 



V 


1 °J 


dw _ 

dt dx ’ 


The initial condition is w(i,0) = f(z) and the boundary conditions are 

u(0,t) = g^t), u(l,t) = g 2 {t). 

A semidiscrete approximation is 

dwj B 


dt 2Ax 


(Wy +I 


(3.49) 


(3.50) 


(3.51) 


The boundary conditions (3.50) are assumed to be homogeneous: 

u 0 = 0, uj = 0. (3.52a, b) 


The NBSs are 


dv o _ Ui — uo dvj _ uj — Uj _ i 

dt Ax ’ dt Ax 


(3.53a, b) 


Gustafsson [3] used this semidiscrete approximation as an example of a problem 
that is Lax-Richtmyer stable for homogeneous boundary data but unstable for highly 
oscillatory non-zero boundary data. We demonstrate that the approximation is GKS 
unstable because of the presence of a a stationary mode. 
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One can show that the NBSs (3.53a,b) are equivalent to the extrapolation formulas 



h(E)u - 1 = 0, 

p{E)u J+ i = 0 

(3.54a, b) 

where 

h{E) = 1 + E 2 , 

p(E) = 1 + E~ 2 . 

(3.55a, b) 

In deriving (3.54) 

we have used (3.52). The 

polynomials h(ic) and p(/c) are 



h(jc) = 1 + K 2 , 

p{k) = 1 + 1/k 2 

(3.56a, b) 

and, therefore, 

h(— 1/k) = h(K.)/t c 2 , 

p(-l/*c) = p(k)k 2 . 

(3.57a, b) 

Futhermore, 

h(— K.) = /i(k), 

p{k) = h(l/*c). 

(3.58a, b) 


As a consequence of (3.57) and (3.58) there is a stationary mode for the finite-domain 
problem. 

There is a close connection between this example and previous example. In fact, 
the example of this section can be written as two uncoupled semidiscrete problems of 
the form given by example 3. The NBS (3.41) is inconsistent but the analogous NBSs 
(3.54) only appear to be inconsistent, i.e., they are actually consistent. 

As in the example of the previous section the eigensolutions of the finite-domain 
problem can be determined analytically. We look for a solution of (3.51) of the form 

wy = c at w y . (3.59) 

The components Wy of the eigenvector w for the finite-domain problem are given by 
Wy = aK?4>+ + 6(-1/k) j ^ + -I- c(-ac) j ^_ + d(l//c) J ^_ (3.60) 


where 



(3.61a, b) 


The constants a,b,c,d are determined by inserting (3.60) into the boundary conditions 
(3.52) and (3.54). One obtains a homogeneous linear system and a nontrivial solution 
exists if and only if 

h{> c) = 0, k 2J = 1. (3.62a, b) 

The polynomial p(/c) does not appear here because of identity (3.58b). 

We first consider h{t c) = 0 where h(/c) is defined by (3.56a). The roots are k = ±i 
and for J even these roots also satisfy (3.62b). For J odd one can show that the roots 
k = ±* lead to trivial functions. From (3.62b) k is a root of unity and hence 


K = t'* llJ , £=l,2,--,2 J. 


(3.63) 
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The associated eigenvalues are 


* ^ — V 

S£ = tsin^— J, <=1,2, •••,2 J (3.64) 

where s = sAi. For J odd, the eigenvalues (3.64) are identical to the eigenvalues of the 
circulant matrix associated with the IVP. Furthermore there are no eigenvalues $ = ±i. 
But for J even there are are two distinct eigenvalues i = ±t. 

One can prove by the energy method that the semidiscrete approximation of this 
section is Lax-Richtmyer stable. But the approximation is GKS unstable because there 
is a stationary mode. The GKS perturbation analysis is identical to the analysis at the 
end of example 3, i.e., an e perturbation in k yields an e 2 perturbation in 5R(s). 

4. CONCLUSIONS 

Stationary modes for semidiscrete approximations are easy to detect because g(/c) and 
q(-l//c) have a common factor. For simple approximations the k roots in class {III) 
defined by (2.1) can be determined analytically and consequently so can the normal 
modes for the finite-domain problem. This includes, of course, stationary modes for 
non dissipative approximations where \k\ = 1 independent of J. Roots in class (7) and 
(II) defined by (1.13) cannot be determined analytically and the corresponding normal 
mode analysis is analytically intractable on the finite domain. 

If there is a stationary mode for a nondissipative approximation on a finite-domain 
problem, then the GKS perturbation test will always indicate the presence of GKS 
generalized eigenvalue. If in the GKS perturbation test the perturbation in k is the 
same order as the perturbation in 3?(s), then the approximation is Lax-Richtmyer 
unstable in the Z <2 norm. On the other hand, if the the perturbation of k results in 
a higher order perturbation in 5R(s), then there is a trivial stationary mode and the 
approximation is Lax-Richtmyer stable. 

If a semidiscrete approximation with a stationary mode is Lax-Richtmyer unstable, 
then the solution exhibits algebraic growth as the mesh is refined. This is a weak type 
of instability and the approximation can be thought of as only marginally unstable 
since the approximation converges in the La norm for smooth initial and boundary 
data on a finite domain. 
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